Cosmological Magnetic Fields from Primordial Helical Seeds 
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Most early Universe scenarios predict negligible magnetic fields on cosmological scales if they are 
unprocessed during subsequent expansion of the Universe. We present a new numerical treatment 
of the evolution of primordial fields and apply it to weakly helical seeds as they occur in certain 
early Universe scenarios. We find that initial helicities not much larger than the baryon to photon 
number can lead to fields of ~ 10" 13 G with coherence scales slightly below a kiloparsec today. 
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INTRODUCTION 

The origin of galactic and large scale extragalactic 
magnetic fields (for which there is no detection yet on 
scales larger than mega-parsecs) is one of the main unre- 
solved problems of astrophysics and cosmology 0]. In 
most scenarios where magnetic fields are produced in 
the early Universe, these seed fields are concentrated 
on scales below the horizon scale where they dissipate 
quickly, and are too small on cosmological scales to have 
any observable effects. However, if pseudoscalar interac- 
tions induce a non- vanishing helicity of these seeds, such 
as in string cosmology |2| or during the electroweak phase 
transition by projection of non-abelian Chern-Simons 
number onto the electromagnetic gauge group ^J, ^j, 
then part of the small scale power can cascade to large 
scales and produce observable effects ]|, |[ [jj ||. In this 
paper we develop a new numerical approach to treat such 
non-linear cascades up to zero redshift and apply it to he- 
lical seed fields produced in the early Universe. 



MHD IN THE EARLY UNIVERSE 

The principal equations for magnetic field B and ve- 
locity field v in the one-fluid approximation of magneto- 
hydrodynamics (MHD) are M 



d t v + (v • V)v = 



d t B = Vx(vxB-i)VxB) 
(V x B) x B 



47T/9 



(1) 



where r\ is the resistivity and p is the fluid density. The 
second equation describes backreaction of the magnetic 
field on the flow. To eliminate it, following Ref. [Q, (To] , 
we write 
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(V x B) x B. 



(2) 



where (3 is the drag coefficient and r = 1/(3 is 
the fluid response time to the Lorentz force f = 
[(V x B) x B]/(47r). The latter can be viewed as the 
time the charged fluid can be accelerated until it inter- 
acts (scatters) with other particles in the background and 
therefore describes damping of the magnetic field modes. 



Again following Ref. (hJ, we express the magnetic 
field in terms of correlation functions My (r, t) — 
(Bi(-x.,t)Bj(y,t)), where r = |x — y|, assuming isotropy 
and homogeneity, 
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where Ml, Mjy, and H are longitudinal, trans- 
verse, and helical magnetic correlation functions, re- 
spectively. Mm = d r (r 2 M£)/(2r) is not indepen- 
dent because of V • B = 0. We define the mag- 
netic field and gauge invariant helicity power spectra 
per logarithmic wavenumber interval b 2 (k) and h{k) by 
E M = (B 2 (r)) /(8tt) = /+ 00 dkb 2 {k) / (Sirk) and H K 



' M 



(B(r)-A(r)) = f Q °° dkh(k)/k, with A the vector po- 
tential, B = V x A. One can show that Ml and H are 
related to these power spectra via 



M L (r) = 



+ 00 



dk ji (kr) 



b 2 (k) 



k kr 

dk ji(kr)h(k) , 



(4) 



where ji(x) — sin(x)/a; 2 — cos(x)/ir. is the first order 
spherical Bessel function. In terms of the usual Fourier 
transforms B(k) = / d 3 r/(27r) 3 / 2 exp(ik • r)B(r) etc., 
b 2 {k) = 47rfc 3 B 2 (fc) and h{k) = k 3 J dfi k B(k) • A(k). 
Eq. (0) also shows that M L (0) = 8ttE m /3, and H M = 
—3 Jq +0 ° drrH(r), and \h(k)\ < b 2 (k)/k implies for all r 

\H{r)\<\M L (r)\/r = H m ^{r). (5) 

Cosmological expansion can be taken into account by 
redefining M L -> M L /T 4 and H -> H/T 5 from now on, 
where T is the cosmological temperature. Assuming for 
now the absence of any external source terms such as 
fluid motions except the one induced by the magnetic 
field, i.e. using Eq. the MHD equations (|l|) reduce 
to 

d t M L = -^d r (r 4 K<9 r M L ) - AaTH 
d t H = \d r [r 4 d r (2 K H + aM L /T)] (6) 



where 
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a = H(0,t), 
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(7) 



and all quantities appearing here are in physical (not 
co-moving) coordinates. The diffusion term n consists 
of a microscopic (rj) and a non-linear drift contribu- 
tion, whereas the a effect is only due to non-linear drift 
here. The source terms will be discussed in the next 
section. If the spatial derivatives of Ml and H fall off 
faster than 1/r for r — > 00, Eq. (g) implies dtHm = 
9 [2kH (0) + aM L (0)] which, together with Eq. (g), shows 
that helicity is conserved in the absence of resistivity. 

Eqs. @ describe small and large scale dynamos of heli- 
cal magnetic fields including damping by Ohmic dissipa- 
tion and " Silk" damping (which is expressed by the red- 
shift dependent relaxation time r) on a unified basis. In 
the early Universe the resistivity can be estimated by r/ ~ 
1/(40ttT) before photon decoupling, T > 0.25 eV @, 

and by the Spitzer resistivity rj ~ 7rm e /2 e 2 /T e 3/2 (where 
m e , e, and T e ~ 10 6 K are electron mass, charge, and 
temperature, assuming full ionization) after recombina- 
tion [|| (the results are insensitive to the latter). Below 
e + e _ annihilation at T ~ 20keV, within the MHD one 
fluid approximation the fluid coupled to the magnetic 
field is well represented by the tightly coupled remaining 
free electrons and protons and r is governed by Thom- 
son scattering of photons off electrons. In this regime we 
use ill 
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(8) 



where the number of free electrons per nucleon X e is ~ 1 
for T > 0.25 eV and ~ 10~ 5 for T < 0.25 eV, and n b h 2 is 
the baryon density in terms of the critical density times 
the Hubble constant in units of 100 km/s/Mpc today. For 
T > 20keV we can approximate the fluid to consist of 
the electromagnetically interacting particles and r is gov- 
erned by neutrino scattering with [I4] 



10 



u /MeV\ 5 8.75 g r 



cm , 



(9) 



and p/T 4 ~ gfir 2 /30, where g r , gf, and g v — 5.25 are the 
statistical weights of all relativistic particles, the particles 
in the fluid and of the neutrinos, respectively. 



HELICAL SEEDS 



Here we consider helical fluid motion, as it can arise 
during cosmological phase transitions (see, for example, 
Ref. H Q for the electroweak phase transition). This 
case has already been treated in Ref. p0| which we adapt 



here to our situation. Since the backreaction of the mag- 
netic field onto the fluid motion has already been taken 
into account by the approximation Eq. (||), the exter- 
nal fluid flow v e is assumed to be uncorrelated with the 
field. It is furthermore assumed that the correlation time 
of the external velocity field is much smaller than the 
time scale of change of the magnetic correlation function, 
(v ei (x, t)w e ,j{y, sj) = Tij(r)5(t - s), where, in analogy to 
Eq- @) r=\x-y\, and 
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(10) 



Assuming for simplicity an incompressible fluid, V • v e = 
0, the additional terms in Eqs. (||) and ([?]) are given by 
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Tl{t) 
- 2C(r) 



(11) 



such that k and a obtain a scale dependent tur- 
bulent diffusion and a effect contributions, respec- 
tively, from the fluid. Here Tl and C are given 
by T L {r) = r corr (r-v e (0)r-v e (r))/r 2 and C(r) = 
T con- ( r • v e(0) x v e (r)} /(2r 2 ). The correlation time r corr 
is either the damping time scale due to interactions with 
the background or, if all components are tightly coupled 
and move as a whole, the age of the Universe t u (T) at 
the relevant epochs. 

The spatial velocity correlations Tl and C can be ex- 
pressed in terms of their power spectra v 2 (k) and c(fc), 
respectively, in complete analogy to Eq. ([|). In general 
they will have the form 

1 



T L (r) = -r COII (v 2 e )(T)f(r) 
\C(r)\ < \T L (r)\/r = C max (r), 



(12) 



where f(r) is a dimensionless function with f(r) — » 1 
for r — * and, typically, a power law fall-off at large 
distances, and the total power (v 2 ) (T) typically peaks at 
a certain temperature T p h, for example, at a primordial 
phase transition, and becomes negligible for T 3> Tph 
and T <§C T ph . 



NUMERICAL SIMULATIONS AND RESULTS 

For any early Universe scenario the initial conditions 
for Ml and H at the temperature where the fields are cre- 
ated should be calculated from the power spectra b 2 (k) 
and h(k), using Eq. @. The magnetic field evolution 
can then be obtained by numerically integrating the non- 
linear partial differential Eqs. (Q) and their extensions 
with helical source terms in co-moving coordinates from 
this initial time up to redshift zero. This is done by 
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employing an alternating implicit method |D| to a one- 
dimensional grid of typically a hundred bins roughly log- 
arithmic in co-moving distance between the inverse of 
today's cosmic microwave background (CMB) tempera- 
ture To and ~ 10 4 Mpc, and using the logarithm of the 
temperature InT as independent variable, adopting the 
standard relations between time and temperature, see 
e.g. Ref. O]. In order to assure that induced velocities 
Eq. (Q) remain non-relativistic, the induced contributions 
to the coefficients k and a, Eq. (M), are limited to the cor- 
responding contributions, Eq. (Thy), of a maximally strong 
external fluid flow during the simulations. 

At a physical length scale r at cosmic time t the ac- 
curacy requirement on the step-size is |H| A In T < 
r 2 / [tMax(K, a)]. For about 10 4 time steps per decade in 
T and for the coefficients given by Eq. M) this is typically 
fulfilled for temperatures up to close to the GUT scale 
and co-moving lengths down to the parsec scale which are 
mostly of interest here. Although this accuracy require- 
ment is not fulfilled at the smallest length scales close 
to the inverse temperature used in the numerical inte- 
gration, the implicit method assures at least convergence 
toward the equilibrium solution at such scales. 

The power spectra b 2 (k) and h(k) can be obtained by 
inverting the transformations of Eq. (0), but a rough 
estimate is given by b 2 (k) ~ M^l/k) and h(k) ~ 
H{\/k)/k 2 . 

In the following we parametrize the magnetic seed field 

by 

M L {r)=N^- 1 (13) 
90 (1 + r/rB) 

where N characterizes the strength relative to thermal 
density, tb is the scale on which it is concentrated, and 
n is the power law index at much larger scales (causally 
produced fields correspond to n > 5). 

To demonstrate the general effect of helicity we start 
with magnetic fields of non-vanishing helicity in the ab- 
sence of source terms. We start at the electroweak scale, 
T = 100 GeV, with a seed field Eq. (|l|) with N = 0.1, 
concentrated at a scale tb = 10~ 4 t M (100 GeV), and 
a power law index n = 3, motivated by a superposi- 
tion of magnetic dipoles, as may be expected for the 
electroweak transition We assume an initial he- 

licity Hm ~ lOOn;,, somewhat larger than the baryon 
number rib, as suggested by models j|, ||. Assum- 
ing the relative helicity, h{r) = H (r) / H max (r) , to be 
roughly independent of r, this corresponds to H(r) ~ 
100/N(n b /n-y)M L (r)/r ~ 5 x lQ- 7 H max (r), where the 
baryon to photon ratio nt/n-y ~ 5 x 10~ 10 . Fig. [I] shows 
results for Ml and the relative helicity. The field at zero 
redshift is decreased by dissipation up to the ~ 0.1 parsec 
scale, whereas inverse cascades enhance the field on scales 
of a few parsecs, reaching ~ 10~ 14 G. For comparison 
Fig. |l| also shows the larger enhancement of Ml for max- 
imal initial helicity (the case discussed in Ref. S) which, 
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FIG. 1: Results in co-moving length scale r for the case 
without external fluid flow. Thick lines show Ml(v) in units 
of T 4 , for initial condition at T = 100 GeV (thick dashed), 
and at zero redshift (thick solid, i.e. Ml = 1 corresponds to 
a field strength of ~ 1.4 x 10~ 6 Gauss today, note that Ml is 
quadratic in B). Thin lines show helicity relative to maximal, 
h(r) = H(r) / H max (r), for initial condition (thin dashed), and 
at zero redshift (thin solid). For comparison, the thin dotted 
line shows the final Ml(t) for maximal initial helicity (not 
shown) . 



however, we consider speculative in the absence of a spe- 
cific model predicting such large helicities. 

It is easy to show that the total helicity Hm which 
is dominated by the peak of h(r) in Fig. [l] is roughly 
conserved, and thus evolution is dominated by magnetic 
back-reaction onto the fluid. Indeed, conservation of Hm 
is usually employed to estimate the field strength via 
B 2 ~ Hm/I c which requires an analytical estimate of 
the coherence scale l c ||. In our numerical approach l c 
comes out without further assumptions as the scale where 
the correlation function cuts off. 

Another interesting situation could be the produc- 
tion of baryon and lepton number comparable to unity, 
rib/riy ~ 1 at T — T n 3> 100 GeV, for instance during a 
phase transition related to new physics, which could give 
rise to maximally helical flows as well. These flows would 
consist of the tightly coupled electroweak plasma and 
could survive as a small perturbation at least down to the 
neutrino decoupling temperature, i.e. i corr ~ t u (T) for 
T > MeV. Their amplitude can be estimated by the dilu- 
tion factor (y 2 ) ~ fafcA^today ~ 10~ 12 due to the neces- 
sary entropy production above the electroweak scale. As- 
suming a causal flow with power on scales not too far be- 
low the horizon scale, we use f(r) = {l+rT/t u (T n )/T n ]~ 5 
with T n = 10 13 GeV for the other parameters in Eq. (O). 
We start with the same magnetic field produced at the 
electroweak transition as above, but this time with van- 
ishing initial helicity, H(r) = 0. Fig. || shows that in 
this case the magnetic field develops helicity and reaches 
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FIG. 2: Results in the presence of a fully helical external fluid 
flow of small amplitude, (vl^ = 10 -12 , at the electroweak 
transition. See Fig. m for line key and text for details. 



where fluid viscosity can not be neglected. The fields we 
obtain are certainly larger than from "astrophysical" seed 
field mechanisms such as the Biermann battery, but are 
also well within the limits from big bang nucleosynthe- 
sis and the CMB (the best of which are ~ 1CP 9 Gauss 
on kpc-Mpc scales, see, e.g., Jl?], O, and from 

gravity wave production (B(r) < 10 1 (r/100pc) -3 for 
n = 3 |p0[). Such fields may also be testable by ultra- 
high energy cosmic ray deflection . The approach pre- 
sented here can also be applied to other magneto-genesis 
scenarios with pseudo-scalar seeds such as in string cos- 
mology Q where coupling to axions may lead to larger 
helicities. 
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values close to 10 _13 G up to about 100 parsecs. The 
coherence scales are also consistent with analytical esti- 
mates 0, |8) , but are considerably larger than in Ref . [|| . 

Our results also demonstrate that the presence of he- 
licity prevents complete dissipation of the fields at small 
scales, resulting in a flat correlation function up to the 
coherence scale. Furthermore, the relative magnetic he- 
licity rises linearly with r and is close to maximal at the 
coherence scale. This could have ramifications for the 
actual detection of helicity, for example, via its effects 
on the CMB [|l6] and could be an important signature of 
physics at or above the electroweak scale. 



CONCLUSIONS 

We used the evolution equations for the two-point cor- 
relation function of helical magnetic fields in MHD ap- 
proximation including magnetic diffusion, fluid viscosity, 
and back-reaction onto the external fluid to evolve weakly 
helical fields produced in the early Universe up to today. 
We find that magnetic fields and/or fluid flows with a 
helicity relative to maximal not much larger than the 
baryon to photon number ~ 10 -9 , as expected during 
the electroweak period, can lead to significant inverse 
cascades. This results in magnetic fields that can be en- 
han ced by several orders of magnitude compared to the 



mer ely red-shifted and frozcn-in initial fields at scales in 



the parsec and kilo-parsec range today. If the seeds are 
roughly thermal in strength and if their power is concen- 
trated on scales not much smaller than the horizon scale 
around the electroweak transition, the coherence length is 
close to a kilo-parsec with field strengths up to 10~ 13 G. 
While this is smaller than the analytical estimates in 
Refs. [Q, it is based on the more realistic assumptions 
of small helicities of order the baryon to photon number 
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